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Abstract 

An effective non-local quantum field theory is constructed, which describes interaction of 
pomerons in the high-coloured QCD. The theory includes both splitting and merging triple 
pomeron vertexes and diagrams with pomeronic loops. The Schwinger-Dyson equations for 
the ’physical’ pomeron are written. Conformal invariance allows to reduce the theory to the 
old- fashioned Gribov pomeron theory with an infinite number of pomerons, one of which is 
supercritical. 

1 Introduction 

The high-energy behaviour in QCD with a large number of colours N c is described by the 
exchange of hard pomerons, which split and merge by a triple pomeron vertex. Exchange of 
coloured objects (single gluons) is damped by factor 1 / . In interaction with heavy nuclei 
the leading contribution comes from diagrams without pomeronic loops (tree diagrams). 
Summation of these is achieved by a closed equation for DIS (the BK equation [1, 2, 3]) or a 
closed pair of equations for nucleus-nucleus scattering [4]. Some estimates of the contribution 
from loops were made in [5, 6]. To take into account the contribution from pomeronic loops 
in a consistent manner one has to consider the effective pomeron field theory introduced in 
[4] as a full-fledged quantum theory. The present study is devoted to this aim. 

Note that pomeronic loops have lately been also actively studied in the framework of the 
colour dipole picture in the so-called JIMWLK approach (see e.g. [7] and references therin). 
There evolution in rapidity of a state considered as a functional of the gluon field is governed 
by a certain Hamiltonian made of the field and functional derivatives in the field. Taking 
loops into account leads to a Hamiltonian containing functional derivatives up to the fourth 
order [8]. In our approach a quantum Hamiltonian can also be introduced. However it 
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contains functional derivatives only up to the second order and so is considerably simpler 
than in the JIMWLK approach. 

In this study we restricted ourselves to a rather formal treatment of the perturbative 
diagrams for the pomeron interaction with or without loops. We construct the relevant 
Schwinger-Dyson equations for the full pomeron Green function and also discuss their con¬ 
formal (Moebius) invariance, which hopefully may simplify their analysis. In fact we are not 
very optimistic about a realistic calculation of the amplitudes with pomeronic loops included. 
This problem presented enormous difficulties already for the much simpler old-fashioned Gri¬ 
bov local supercritical pomeron model, so that even its internal inconsistency was claimed [9] 
(see also [10] for a discussion of this inconsistency). Possibly the more complicated structure 
of the BFKL pomeron and its interaction may overcome these old troubles. However in the 
present study this problem is not touched but left for future investigations. 

2 Effective action and diagrams 

2.1 Effective action 

Our main tool will be the non-forward BFKL Green function as a function of gluon coordi¬ 
nates and rapidities, which satisfies the equation 

(| + H)g{y- y'-,ri,r 2 ;r' 1 ,r' 2 ) = S(y - y')^ 2 ^ 2 5 2 {r iv )8 2 {r 2 2') % (1) 

where r\y = r\ — r\ etc and H is the BFKL Hamiltonian [11]: 

H = ^(ln p\ + In p\ + 4 In ^2 'Pi + 4 lnr i2 ' v\ ~ ^l 1 )) ( 2 ) 

z \ pf p 2 ' 

and d = a s N c /ir. To economize on notations, in the following we shall denote as 2 the point 
in the space formed by rapidity y and two transverse vectors r\ , r 2 : 

2 = {y,ri,r 2 } = {y,p}. 

In this notation we write g(y — y'\ rq, r 2 ‘, r(, r' 2 ) = g(z,z'). The Green function g(z,z') is 
invariant under conformal (Moebius) transformations of coordinates.. 

The inverse function g~ 1 (z, z') defined by 


dz"g 1 (z,z")g(z",z') = 5{z - z') 


(3) 
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where dz = dyd 2 r\d 2 r2 = dydp and 5 (z) = S(y)5 2 (ri)5 2 (r2) = 5 (y) 5 p, is however not confor¬ 
mal invariant, since the measure is not conformal invariant. It is not difficult to construct a 
conformal invariant inverse g~^ v using the invariant measure: 


Then we can rewrite (3) as 


dr = 


d 2 r\d?T2 


12 


(4) 


J dy"dr"r\ 2 g 1 (z, 



= rf 2 S ( z ~ 


(5) 


This shows that the conformal invariant function is 


9inv 


(z, z') = r\ 2 g ( z,z')r' 12 


( 6 ) 


Now we pass to constructing the effective non-local quantum field theory, which is to 
describe propagation of free pomerons and their triple interaction. It has to generate all 
diagrams built from the pomeron propagator, which is the Green function g, and triple 
pomeron vertex introduced in [12, 13, 14]. Let (j> B (z) and 4>a(z) be two bilocal fields for the 
incoming and outgoing pomerons. To reproduce the correct propagator the free action So 
has to be 

So = J dzdz'(j)A{z)g^ 1 (z,z')(j)B{z). (7) 

We introduce interaction with an external sources as 

S E = ~ j dz(f>A(z)J B (z) +(a<t+B^. ( 8 ) 

Here the sources for the projectile A and target B are different from zero at rapidities y = Y 
and y = 0 respectively: 

Ja(z) = JA(p)d{y ~ Y), Jb(z) = J B (p)S(y). ( 9 ) 


Finally the interaction to reproduce the 3P vertex at large N c has to be taken in the form: 

s I = 2 -^f + (A»B), (10) 

71 J r 12 ? 23 ? 31 V 7 

where z\ = {y,r2,rs}, z 2 = {y,r 3 ,ri}, Z3 = {y, r\,r2} and the conformal invariant operator 
L 12 is 


L12 = rf 2 V 2 V 2 . 


(11) 
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Note that the form (10) assumes the fields to be symmetric in the two space points rq and 
r 2 - In fact the symmetry properties of the fields are determined by symmetry properties of 
the external sources. We assume them to be symmetric under the interchange rq <-► r 2 - 

This action leads to BFKL pomeron diagrams with the standard 3P interaction in the 
presence of an external field. Note that the signs of different building blocks of the diagrams is 
somewhat different from the standard ones: the propagator g and external sources enter with 
a minus sign as a consequence of the choice of signs in the action. This latter corresponds 
to the desire to make the interaction real and not pure imaginary as in the original Gribov 
reggeon field theory. 

In absence of the external sources this action is explicitly conformal invariant provided 
(f>B and 4>a are invariant. Indeed the free part can be rewritten as 

So = j dydy'dTdT , (j) A {z)rl 2 g~ 1 {z,z')r' 12 i 4)B{z) (12) 

and the interaction part as 

Si = J d-ydy'dy"drdr'dr"(\> B [z')4>b{z”)^{z!, z"\z)Li 2 ^a{z) + (A B^j , (13) 

where 7 (z 1 , z"\z) is the bare interaction vertex for the incoming pomeron at z and two outgoing 
pomerons at z' and z": 

l{z', z"\z) = S(y - y')5(y - y / ')^ 2 (^i 2 ')^ 2 ( r i , 2 ")^ 2 ( r i" 2 )ri 2 ri, 2 ,ri,/ 2 //. (14) 

This is a conformal invariant function. Since both r\ 2 g~ l {z,z')r' 12 A and 7 ( 2 /, z" | z) are con¬ 
formal invariant, so is the action Sq + Sj. Of course in physically relevant cases conformal 
invariance is always broken by the external sources, which also introduce a mass scale into 
the theory. As a result contribution from any Feynman diagram without external sources 
is conformal invariant, which becomes explicit if one uses the invariant integration measure 
dz{ nv = dydr and invariant interaction vertex 7 given by Eq. (14). 

Finally note that from (1) it follows that as an operator in the 2 -space 

„-'(*,/) = v? (A+ F). ( i5) 

Here we use the notation 

v 2 = V?V|. (16) 

Note that operator H is symmetric: 

V 2 H = | (p\pl In [plpl) + pi In r\ 2 p\ + p\ In r\ 2 p\ - 4^(1 )p 2 lP £). (17) 
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As a result So can be written directly in terms of the BFKL Hamiltonian: 

So = J dydp(j) A (y,p)V 2 p + <j) B (y,p) = j dydp<j> B {y, p)V 2 p + 4 >a{v,p)- 

V V ( 18 ) 
In this form the symmetry between the projectile and target is made explicit: it has to be 
accompanied by changing y —> —y. 

Using a representation for H ( [15]) 

Hf(n,r 2 ) = f d 2 r3r 2 12 (/(n, r 2 ) - f(r i, r 3 ) - /(r 2 , r 3 )) (19) 

J r 13 ? 23 

we can rewrite the free part of the action in a more explicit form 

S o = J dyd 2 rid 2 r2(t>A(y,ri,r 2 )'V 2 p ^-^-(t) B (y,ri,r2 ) 

_a I U B (y,r 1 ,r 2 ) _ (/> B (y,ri,r 3 ) ~ 0 s( 2 /,r 2 ,r 3 ))}. ( 20 ) 

Z7T J r 13 r 23 

2.2 Diagrams and their order of magnitude 

The quantum field theory described by the action S = Sq + Si + Se allows to construct 
the perturbation theory by the standard technique, expressing the amplitude as a sum of 
Feynman diagrams. It is instructive to see orders of magnitude of different contributions in 
terms of two independent small parameters of the theory, d = a s N c /ir and 1/-/V C . Obviously 
each triple interaction contributes a 2 /N c . The final estimate depends on the magnitude of 
the external source. If we treat it perturbatively to be consistent with the whole approach 
then it should correspond to the quark-antiquark loop, from which we extract its order, which 
is a. As a result, the diagram with external lines (sources) and L loops has the order 

< 2i > 

As one sees, for a given number of external sources, the dominant contribution comes from 
the tree diagrams, each loop introducing a small factor a^/N 2 . 

However this pure counting does not take into account the growth of the pomeron prop¬ 
agator at large y as exp Ay where A is the BFKL intercept, nor the enhancement related 
to the nuclear sources with large atomic numbers. For the AA collision amplitude with the 
overall rapidity difference Y this changes (21) to 

(a^/ay(A) 1 2L+ ‘ B ~\^y. ( 22 ) 
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where n is the maximal number of the exchanged pomerons at a given rapidity. This number 
depends on the topology of the diagram and generally grows with Ie and L. 

So for 

cL 4 1/3 ~ 1 ( 23 ) 


one has to sum all tree diagrams and for rapidities Y such that 


a2 

T-e ~ 1 

N r 


(24) 


one has to sum also at least some of the loops. For onium-onium scattering ( Ie = 2) a 
convenient method to sum the leading contribution is to join two sums of tree diagrams 
starting from the projectile and target at mid-rapidity [5]. 

However all these perturbative estimates are to be taken with caution. The full pomeron 
Green function (sum of the so-called enhanced diagrams) may have an asymptotic behaviour 
at large y very different from the bare one. In fact there is every reason to believe that the 
former will grow at most as a power of y. not as an exponential. Then all the above estimates 
will have to be reconsidered. Unfortunately at present a reliable summation of all enhanced 
diagrams cannot be performed even for the old Gribov local supercritical pomeron model. 


3 The pomeron Hamiltonian and operators 


With the action fixed, one can easily construct a Hamiltonian formulation for the state 
evolution. Since the Lagrangian is of the first order in derivatives in rapidity the Hamiltonian 
is just the action without the derivative terms with a minus sign and integration over y 
dropped: 

n = - J dpd> A (p)v 2 ,n<pB(p) 

_ 2 otK + . ^ s } _ 

^ J r 12 r 23 r il 1 V 7J 

where P 23 = {r 2 ,r' 3 } etc. To pass to the quantum theory one has to consider 1 j>A,B{p) as 

operators. Their commutation relation can be easily established from the form of the free 
Green function 

g(y - y'\p,f>) = - < T{^>B{y,p) 4 >A{y\p')} > • (26) 

From (26) we conclude 


r\ 

(dy +H ) K T ^ B ^ y ' p ^ A ^ y ^ p '^ >= ~ y')\ ( t ) B(y,p),(l>A{y,p% 


(27) 
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where we have used the equation of motion for <f>B- Comparison with (1) gives 

[4>B(y,p),4>A(y,p')] = -V" 2 <S(/9- p'). (28) 

Thus in a representation in which 4>b(p) is diagonal and the state vector is a functional 
^>{4 >b(p)} the field 4 >a is essentially a functional derivative 

Mp) = v p 2 5(j)B ( p y (29) 

The state with a given field 4>a(p) will be represented by a an exponential 

*<Mp)({<M) = e /dp<Mp)V?<Mp) - (30) 


The state vector will satisfy the evolution equation 


dy 


HV, 


(31) 


where in Hb, given by (25), one has to substitute the field 4>a by functional derivatives. In 
this substitution, as always, the order of the operators is actually undetermined. If we put 
all the derivatives to the right (’normal ordering’) then explicitly 


2vr . 


d 2 rid 2 r2d 2 r3rl 2 


| <t>B{pv2.) ~ <Pb(pi 3 ) - (Pb(P23) 


d(t>B{pv 2 ) 


+ Ll2fe<m) ( V;2 5tek))( V;;2 5teI^))l}' (32) 

One can pass from this ’target’ representation in which the field 4>b is diagonal to the 
’projectile’ representation in which it is 4 >a which is diagonal and 4 >b is represented by a 
functional derivative 

= “ v ' 2 ^fw (33) 

In this representation the Hamiltonian will be obtained from (32) by interchanging A and 
B and changing signs of derivatives. The state vector will be obtained by a quasi-Fourier 
transformation using (30). 

One can construct a formulation in which the symmetry between target and projectile is 
more explicit. To do this one can pass to slightly different field variables for which the BFKL 
Hamiltonian becomes symmetric. The form of these new variables is clearly seen from the 
commutation relation (28) 


<t>{y,P) = y/^AM = T 4>A{yiP), 



(34) 


4> ] (y,p) = \f^i ( t ) B{y,p) = T4> B (y, p). 

For them the equal rapidity commutation relation takes the form 

[4>{y,p),4> ] {y,p') = $(p - p')- (35) 

One can also assume that the scalar product of state vectors is chosen to make (j) and <jy 
Hermitian conjugate to each other. In the representation in which say (jy is diagonal with 
complex eigenvalues a we take (up to a normalization factor) 

<^ 1 |^ 2 >= j'DaDa*^ 1 (a*)^ 2 {a)e~ f dpa ^ p)a(p) . (36) 

Then indeed 

< Ti^(, 0 ) 1^2 >= J DaDa* ^ 1 (a*)^ 2 (a)a(p)e~ -f dp ' a * (p ' Mp,) 

= J DaDa* 'F 1 (a*)'F 2 (a)( - ^-y)e“ / d P' a *WW) 

= [ DaDa* ^ 2 (a)e~ I dp ' a * (p ' )a{p ' \ 8 . . fri (a*) =< <j>(p)V i |* 2 > • 

J oa*(p) 

Thus the two quantized fields cj)(y, p) and (f)'(y, p) acquire the standard meaning of annihilation 
and creation operators for a pomeron at rapidity y and space points p = {ri,r 2 }. 

In terms of these new field variables the free action takes the form 

So = J dydp 4 >(p)T (J?- + Hj T~ l 4 >\p) = J dydp(j>{p) (J^ + H^j ${p ), (37) 

where a new Hamiltonian for the pomeron is 

H = THT~ l = ^(ln p\ + lnp 2 + lnr i 2 \/^ + \f^ lnr i 2 \/^t ~ 4 V ; ( 1 )) • (38) 

2 v V Pi V P '2 V P ‘2 V Pi 1 

It has obviously the same eigenvalues but is Hermitian (and real). Using its hermiticity we 
can revert the order of operators in (37) and write So in the ’normal order’ form 

So = J d/ydfxpHp) ^--^- + 5^ (j>{p). (39) 

This form explicitly shows the symmetry between target and projectile, which is quite 
similar to the usual time reversal: one has to change <f> «-► (fy, y —■► —y and revert the order 
of all operators. 

In terms of new field operators the external part of the action aquires the form 

Se = - J dz(4)(z)T~ l J B (z) + <$ Ja^z)^ = - J dz(<j>(z)j\z) + h.cj, 


(40) 
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where 

J(z) = T~ 1 J a {z), jftz) = T~ 1 J b {z). (41) 

The interaction part becomes rather complicated, involving several operators T or their 


inverse: 


s = 2 a%N c 


7T 


dy d r ] d 2 2< i r3 ( T V t (y.P 23 )-r Vt(y>Pi3) • r\ 2 Tft(y,p V2 ) + fr.c.). (42) 

rfn'rfy.r... V / 


12' 23' 31 

With the physical meaning of operators (f> and ft well established and indeed standard, the 
analysis of evolution becomes trivial. Let us follow it for free pomerons. Then their number 
is conserved and actually the only connected diagram corresponds to a single pomeron. Such 
a state is to be constructed as a superposition of single pomerons at different positions p\ 


®{y)= / dpf(y,p)ft(p)A> o, 


(43) 


where Tq is the vacuum state which obeys 


0(p)T o = 0 (44) 

and is normalized to unity, (we assume the Shroedinger-like picture with operators ft and ift 
at fixed rapidity). At the initial rapidity y = 0 the pomeron wave function is determined by 
the external current: 

T(0) = I dpJ*(j>)<ft(p)* 0 , (45) 

where we recall that J is the spatial part of J. This state evolves to the final rapidity Y at 
which we are interested in the amplitude A to pass to a specific final state determined by 
the current at y = Y: 

47 = j dpJ(p)(ft(p)^ 0 . (46) 

One has 

A f i =< 4'/|4'(Y) >= J dpdp J(p)f (y, p') < ^ 0 \ft{p)cft{p')\^ 0 >= J dpJ(p)f(y,p ), (47) 
where we used (35) and (44). 

The law which governs evolution of the wave function f(y,p ) follows from the general 
Schroedinger equation (31) and the form of the Hamiltonian 7i. The free part of the latter 
in terms of new operators has the standard form 

H 0 = - J dpcft { 


z)Hft(z) 


(48) 
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so that from (43) one immediately finds the equation 

^ (49) 

dy 

with a formal solution 

f{z ) = e~ Hy f(0 ) = J dp'g(y, p- 0 , p’)f{0 ). (50) 

Here g is the Green function for the operator d/dy + H which can be written as an operator 
in the coordinate space 

g(y) = 0(y)e~ Hy . (51) 

Using this we obtain for the amplitude 

A fi = J dpdp' J(p)g(y, p;0, p')fi (p 1 ). (52) 

Returning to the initial external sources and the Green function we reproduce the standard 
result 

Afi = J dpdp'J A (p)g(y, p; 0, p')J B (p). (53) 

Indeed we have 

s=(|; + U"= r (| + ff )" r_1 = r ^ ( 54 > 

Putting this into (52) gives (53). 

4 The Schwinger-Dyson equations for the pomeron Green 
function 

4.1 The pomeron self-mass 

The pomeron self-mass operator starts and finishes with the 3-pomeron vertex, which contains 
operator L acting on the incoming and outgoing pomeron propagator. As a result, the Dyson 
equation for the full Pomeron Green function G(z, z') takes the form 

G(z, z') = g{z\z') - J dz inv dz' inv g(z , z)LE(z, z')LG(z ', z'), (55) 

where both operators L act on pomeron propagators (the minus sign in front of the second 
term is due to the propagator being in fact equal to —g). Applying to this equation operators 
L from the left and from the right we find 

G(z, z') = g(z , z') - J dz inv dz' inv g(z, z)Y,(z, z')G(z ', z'), 


(56) 
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where we define 

G = LGL, g = LgL. 

Eq. (57) can be rewritten in an obvious operatorial form as 


G = g - ~gYG, 


(57) 


(58) 


which is the standard form for the Dyson equation, except for the sign. 

Note that the self-mass £ entering this equation is a conformal invariant function, due to 
conformal invariance of both the pomeron propagator and 3-pomeron vertex. As a result the 
full Green function is also conformal invariant. 

To pass to pomeron ’energies’ lo we have to understand how they are related in the 
3-pomeron vertex. We standardly present 

ra-\-io o rf, , 

G(V) = / (59) 

J a—ioo i 

with the inverse transform 

/ +oo 

dye~^G(y). (60) 

-OO 

In the lowest non-trivial order, suppressing the integration over coordinates and coupling at 
the vertexes, we have 


G(Y) = / dydy'dy"S(y + y' + y" - Y)g(y)g 1 {y'')g 2 {y")g(y') 


(61) 


where g\ and 52 are the propagators of the two intermediate pomerons. For G(ui) we find 
G{u)= f dydy'dy l 'e~^ y+y Gy , ^g{y)g l {y'')g 2 {y l ')g{y') = g(u) [ dye~“ y g\(y)g 2 (y)g{u). (62) 


The integral over y can be written in the form 

J J a-\ — too 1 J ao—100 Ant 


(63) 


>a 1—200 

We can always choose u> on the line 


Re (lo — ui\ — L 02 ) = 0. 


Then integration over y will give 


27n5(Im (lo — oj\ — 002))- 
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Subsequent integration over, say, ui 2 will finally give 

r rai+ioo 

dye wy gi(y)g 2 (y) = — gi(ui)g 2 (u - u 2 ), (64) 

J J a\— ioo 

which result should be analytically continued to arbitrary oj. This means that in terms of 
energies the 3-pomeron vertex formally contains 

2m5{u — ui\ — u) 2 ). (65) 


that is, energies are conserved at the vertex. 

Using this result, again in the lowest order, we find, explicitly showing the coordinates: 
(1,2) = {ri,r 2 } 


S(, 0) (1,211', 2') 


8a gN% r du 1 

7 r 2 J 27 n 


ri 2 d 2 r 3 r' 12 d 2 r ' 3 
r 2 2 / 2 , 2 

'31' 21 r 31 r 21 


<7 aJ1 (l,3|T,3 / ) 5 „(2,3|2 / ,3') 


( 66 ) 


To pass to the full self-mass we have to substitute the full pomeron Green functions G 
for the propagators g and change one of the 3-Pomeron vertices into the full one T. In this 
way we obtain 


£ w (l, 2 |l', 2 ') 



du 1 
27 n 


^|^dr(l",3")dr(2",4") 

r 31 r 21 


G wl (l,3|l , ',3 ,, )G a; _ a 


(2,3|2",4")T^ 1 (1",2",3",4"|1',2'). 


(67) 


In the lowest order the vertex is given by the spatial part of 7 , Eq. (14): 


r(°)(l", 2", 3", 4"|l', 2') = 5 2 {l" - l')<5 2 (2" - 2')«5 2 (3" - 4")r 2 , 2 ,r 2 „ 3 »r 2 2 „ 4 „. 


( 68 ) 


4.2 The vertex equation 

As in any quantum field theory with a triple interaction the full vertex is determined by an 
infinite sequence of skeleton diagrams, containing the vertex itself and full Green functions. 
In the lowest order one has the so-called ’3-gamma’ equation. Its explicit form is rather 
complicated due to large number of variables. From the corresponding Feynman diagram we 
find in terms of gluon coordinates 

r^(T, 2', 3', 4'|1, 2) = r (0 ) (T, 2', 3', 4'|1, 2) 

_8a^Nl r dui f dr ^y'^ 2")d T (3", 4 // )d r (3" / , 4 / ")c?r(l, 2)dr(V, 2')dT(l'", 2'") 
it- J 2 tti J 

L( 1 ', 3')r^ M (T, 3'|1", 2", 3", 4")G W1 (1", 2"|T", 2'") 
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G W '- W1 (3", 4"|3"', 4', 3 W , 4"'|1, 2) 

L(l, 2)G W _ W1 (1, 2|1', 2')r WiW1 (l"', 2'", 1', 2'|1, 2). (69) 

In the shorthand notation p = { 77 , r 2 } and correspondingly 

6 ^( 1 , 2 ) = G UJ (p 1 ,p 2 ) = G(r ( ^ ,r ( 2 1] r[ 2) r ( 2 2) ) 


the equation for the vertex part can be rewritten in a more compact form 

r W)W '(l , ,2 / |l) = r (0) (l , ,2'|l) - J^LJ dr(l // )dr(2 // )dr(3' / )dr(4 / ')dr(5 // )dr(6 // ) 

L(l0r^ jWl (l / |l / ',2'')G wl (l"|6 / ')G a; ,_ a , 1 (2''|3' / ) 
^.- W1 ^(2^3^4'0L(4 / 0G^ 1 (4' / |5 / 0^ w ^ 1 (5 / ^6 / '|l). (70) 


5 Conformal (Moebius) invariance 

5.1 Conformal basis 


The basic building blocks of the pomeron perturbation theory: the propagator g , bare 3- 
pomeron vertex 7 and integration volume dr, - are conformal invariant. As a result also 
the full Green function G, self-mass £ and vertex T are conformal invariant functions of 
their arguments. It seems profitable to use this property to simplify the Schwinger-Dyson 
equations. To do this we have to study a general form for the conformal invariant functions 
for the transitions pomeron—» pomeron (’two-point functions’) and pomeron—» 2 pomerons 
(’three-point functions’). This can be achieved by expanding these functions in the conformal 
basis formed by functions (in complex notation) [16] 

1+ra 


E,(p) = E fl (r 1 ,r 2 ) = 2 + " (^V) 

\rwr 20 J Wo r 20 / 




(71) 


where p = {n,is,ro} = {h, ro} with n integer, v real and two-dimensional transverse ro 
enumerate the basis. 

Functions E fl {p) are the proper functions of the operator L 

47r 8 


LE fl (p) = l lt E lt (p), 1^ = 




(72) 


where 


®n,v — O'ft — 


7T 4 1 

2 v 2 + n 2 /4 


(73) 
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They form a complete system: 


ri 2 6(p-p') = -J2 E M E Z(p'), 


1 


where we use a notation 

OO r 

E = E ^ I A,, 

/x n =—oo 

and satisfy the orthogonality relation 


drE lx (p)E*,(p') 


= a n ^5 nn >8(v-is , )5 2 (roo>) + b n „8 n - n '6(v + v , )\roo'\ 2 ilv (^-) n 


00 ' 


The coefficients b nu are given by the formula 


b n , , = 7T 3 


■)4x^ 


T(-?T + (1 + \n\)/2)T(iv + |n|/2) 
—+ |n|/2 T(zi/ + (1 + |n|)/2)r (—iu + |n|/2) ’ 


Functions E, t are not linearly independent. In fact 


/ r * \ n 
I-2+4*1/ / 1 00' 


E-n,-v,rM = — / Iroo'T™ (p). 

flni/1 \r 0 0'/ 0 


Presenting a function f(p) as 


/(/?) = J dp'8{p- p')f(p') 


(74) 


(75) 


(76) 


(77) 


(78) 


(79) 


and using completeness (74) one gets 

f(p) = \! dr' ^2 Efi(p)E*(p')f (p') = We^p)?,, (80) 

where 

U = J drE;(p)f(p). (81) 

This gives the standard expansion in the whole overcomplete basis. 

However it seems possible to limit oneself to an independent part of this basis. One 
possibility is to take a restricted basis with v > 0(u < 0). We denote this restriction by 
/j > Q(p < 0). In fact we may split the integration over v in (80) into two parts: 

f(p) = Je E ^(p)fn +1J2 E v(p)f^ 


/t>0 


p ,<0 


(82) 
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and, say, in the second term express the eigenfunctions E^ with p < 0 via those with p > 0 
using (78) to obtain for the second term 


1 


/j,<0 


b—n— 


&—n—b 


/<! 2 >> m ,|- 2 - 4 “- ()^7) 


(83) 


Changing summation and integration variables n —> —n and v —* — v we get for this term 

\ Y f-n,-u,ro— ! d 2 r' 0 \rw\- 2+Alv \ E , (p). (84) 

2 il>0 °w J V r oo'/ 0 


Interchanging also ro and r' 0 integrations we finally find for it 

1 
2 


YMp)Ui 

p> 0 


where 


/, = — / d 2 >-;/-„,-,, r ;l>' 0 ' 0 |- 2+4 “' (^V 

a nv J 0 Vr 0 /o/ 


(85) 


( 86 ) 


Summing this with the first term in (82) we get the desired expansion in states with p > 0 

m = \ E Mp)(u + h) = E (87) 

ti> 0 /i> 0 

Integrating this with E* >0 (p) and using (76) we find 

= f drE*(p)f (p). (88) 

Putting this into (87) we find 

f{p) = Y E,{p) [ dr'E*{p')f(p'), (89) 

M>0 J 

which means that one also has a completeness relation for half of the basis with v > 0 : 

ri 2 6(p-p') = Y E »(p)K(p')- (90) 

/i>0 

Obviously the same property is valid for the second half of the basis with v < 0. 


5.2 Two-point functions 

We have to deal with a conformal invariant function T(l|l'), where, as introduced in the 
previous sections, the arguments refer to pairs of the pomeron coordinates. Using half of the 
conformal basis we present A as 

T(i|io= y ^(i)^'(iK' = E E ^(i^'dOv. 

H,H r >0 fJ.>0 i_i '<0 


(91) 
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where we use 

E n ,v,r 0 (p) = E- n -v,r 0 {p)- ( 92 ) 

Our aim is to see which properties A^y should have for the function A(l,1') to be conformal 
invariant. Since E n ^ ro (n,r 2 ) = -En,iv 0 +a(?’i + a, r 2 + a) the translational invariance requires 
A^i to depend only on the difference ro — r' 0 = roo'- Under inversion the expression to be 
summed over n > 0 and [i! < 0 changes as follows 

•VO-coO - . (a. /.) • A m , (^-) . (93) 

Here and in the following (a. f.) means ’antiholomorhic factor’, that is the complex conjugate 
of the preceding factor. Invariance under inversion requires that 

A m , (^) = rJ+"- 2i ‘V; 1+ "'- 2i,/ ■ (a. /.) ■ A m ,(r m ). (94) 

However the left-hand side only depends on the product ror' 0 , so in the right-hand side we 
are obliged to have either n = n' and v = v' or A MM '(roo') oc 5(roo')- Since v and v' have 
opposite signs, the first alternative cannot be realized. With A AlM '(roo') oc <5(roo')- we find 
that n + n' = 0 and u + u' = 0 so that 

A^n' — Sn,—n'$(y A V )fi('roo l ')Q"n.i/A n i/ = 6^A^, ( 95 ) 

where fi = g(n —> — n, —> — v) and we defined = A nu with factor a nv separated for 
convenience. As a result the double sum in (91) transforms into a single one 

A(l|l') = £ E„(l)E p (l')A nu = J2 E,(l)E;(l')A n „ (96) 

fi> 0 /i>0 

where ft = fi(n —> —n, /z —> —v). A similar form with summation over fx < 0 can be obtained 
in the same manner. Representation (96) or a similar one with a sum over g < 0 are valid 
for any conformal invariant two-point function. Note that taking an average of the sums over 
H > 0 and fx < 0 one obtains a similar representation in terms of the whole overcomplete 
basis, which is standardly used for the BFKL Green function g(z, z') 

Now suppose we have conformal invariant functions A?(l11') and (7(1, 1 ; ) and form a 
conformal invariant integral 

A{ 1|1') = I dr ,, R(l|l")<7(l"|l / ). 


(97) 
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Each of the three functions A, B and C has the representation (96) with conformal coefficients 
Afj, Bfj, and C^. Doing the integration with the help of orthonormalization properties of the 
basis functions E^ with /i>0we find 

^(1|10 = E E^E^l^B.C, = E (98) 

H>0 n>0 

which means that in the conformal representation (96) 


Afx — BfjC'ij. 


(99) 


Applying this result to the Dyson equation (58) in the u representation we immediately 
find 


Gojfj, — Quji . 


( 100 ) 


So both the energy oj and conformal quantum numbers g of the pomeron are conserved in 
the interaction and its full Green function in the conformal basis is trivially expressed via its 
self-mass: 


1 


GuJIX 1 /g uu + Sa 


1 


01 1 /g uu + HVu 


( 101 ) 


yujf± \ i v 

In Eq.(lOl) the free pomeron Green function g in the conformal basis is given by 


2 1 
dun — — ' 


( 102 ) 


5.3 3-point functions 

For a 3-point function the expansion similar to (96) reads 


r (l|2,3) = E ( 103 ) 

/U>M2M3>0 

where the intermediate c.m. coordinates are R ,\, R 2 and R 3 . Conformal invariance allows to 
determine the dependence on them of F /tl /i2/J3 . Translational invariance requires r /(1 /J2/J3 to 
depend only on differences Rik and from the scale invariance such a dependence should be a 
power one. So we seek 


= Ri?R%?RT • («• /•) • r ni(/1 | n2i/2n3 , 3 . (104) 

After inversion we find an extra factor in the sum (102) 

^,-2—«12—031+1—ni+2ji/i ^ a j j _ g-iniir 


(105) 
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times two similar factors which are obtained from (104) by cyclic permutations of 123 and 
conjugation. Invariance requires that each power is zero and that the sum m + n 2 + n 3 be 
even. We get six equations to determine the a’s. Their solution gives 

<212 = “ + 7j(n 2 - n\ - n 3 ) + i{v 1 - u 2 + v 3), 

«23 = + ^(ni + n 2 + n 3 ) - i(iq + u 2 + v 3 ), 

1 1 

« 3 i = -- + ~{n 3 - m - n 2 ) + i{u 1 + ^2 - ^ 3 ) (106) 

and similar expressions for the powers in the antiholomorhic factor d’s with rq —> —rq. 

For the free vertex T^ given by Eq. (68) we have 

= ( 107 ) 


where again /./ = n{n —► —n, zz —> — z/) and the vertex Ij il/t2/ , 3 has been introduced by 
Korchemsky [17]: 




M1/Z2M3 


d 2 r\d 2 r 2 d 2 r 3 
r 2 2 2 

' 12 ' 23'31 


E l _ ll (ri,r 2 )E IX 2 (r 2 ,r 3 )E l _ l 3 (r 3 ,r 1 ) = V.{h u h 2 ,h 3 ) J J 

i<j 


-A i 


(108) 


with A 12 = /ii + /i 2 — h 3 etc. 


6 S and F in the conformal basis 

To illustrate the simplifications introduced by transition to the conformal basis in this section 
we study the pomeron self-mass and triple interaction vertex in this basis 

In a shorthand notation for the gluon coordinates 1 = (ri,ri) the pomeron self mass is 
given by 

S w (l|l') = J^J dr(2)dr(3)dr(2')dr(3') 

T (0) (l|2,3)G W1 (2|2 , )G tt ,_ tJl (3|3')r W)Wl (2', 3'|l'). (109) 

The actual number of integrations is in fact smaller due to (5-functions in the conformal 
vertex T^ 0 ) defined by Eq. ( 68 ). We expand both the Green function and the vertexes in 
the conformal basis with fi > 0 (see 96) and (103). In the following we also suppress the 
w-dependence to economize on subindexes. We have 

G(2\2 i ) = Y / G^(2)E;(2'), 

ii> 0 


(110) 
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G(3|3') = £ GsE^EZ'tf), 
n’>o 

r (0) (l|2,3) = £ C ( , ira E„(l)£,*, 1 (2)E; i (3), 

Ml,M2,M3>0 


( 111 ) 

( 112 ) 

(113) 


r (2 , ,3'|i / ) = £ r #4 ^ lMi ^,(i')^(2')^(3'). 

Ml 5M2 5 ^ 3 ^*^ 

Integrations over the gluon coordinates are done with the help of (76) and give the product 
of (5-functions: 


a M a M' ^MM2 ^MMj/ ^m'M3 • 


So we find 


E(i, 2 |i', 2 ') = 8 “ iJV? ' ^ 


7T 


E # ^^^^^( 1 )^( 1 ') 


MUMi.M.M 


= E ^^(i)^^ 1 ') 

Ml ,Mi 


where 


MlMi 


8 a^A(? r dm 1 (o) 


7T 


2jj-^ / -< /ail/Xj/i'^M^M^MiMlMi’ 
M,M' 


(114) 


(115) 


and the suppressed dependence on cm is obvious from its conservation at the vertexes. Thus 
we have found for E(l, 2|1', 2') an expansion in the conformal basis (91). From its conformal 
invariance it follows that 

S /UMi = WEmw ( 116 ) 

where E^ is the desired pomeron self mass in the conformal representation. It can be found 
from (115) after summation over /j,\ This gives 

8 a^N^ f du>i 


v 1 _ w s c 

— 2 
7T Z 


L \' p(o) n r 1 v 
; 1 mImi,m2 mi M2 mi, 


27T* 


M2 | M3 ‘ 


(117) 

Ml.M2,M3 

The sum over //.*, i = 1,2,3 includes integrations over three c.m coordinates i?j on which 
only the vertexes depend. We get an integral depending on the four conformal weights 

Ih\h lM \h 3 = f d 2 R 1( l 2 R 2 d 2 R 3 R^R^R^R^R^R^ • (a./.). (118) 


Here Rq is arbitrary since the integral is independent of it. The powers are given by 

a 0 i = ~ + ^(ni - n 2 - n) + i(v 2 ~v\~ v), 

012 = ~ + ^(ni + n 2 + n) - i{u 2 + f\ + v), 



20 


a 21 = -- - -(ni + n 2 + ra 3 ) + i(i ^2 + vi + v 3 ), 

«3i = + ^(” 2 - n i + n .s) + i(^i - ^2 - ^ 3 ) 

and «02 and 032 are obtained from «oi and 031 by interchange l<->2. The integral (118) is 
convergent both in the infrared and ultraviolet. However its calculation does not look simple. 
Once this integral is known, the self-mass in the conformal basis is given by a sum over three 
conformal weights: 


So/^N^ f doj\ 

^2 J 2lTi ^ ^ (kh)^/l2 ^1 ,fe2|/l3 (^’ 


hi,h, 2 ,h 3 


Here we made explicit the uj dependence introducing it in the arguments; H is the Korchemski 
vertex (107), ^h lt h 2 \h 3 is defined by Eq. (104) and J2h is given by (74) without integration 
over vq. 

In the same manner one can obtain expressions for the vertex part T. We shall limit 
ourselves with the ’3-gamnra’ approximation, Eq. (79). As before we expand the vertex parts 
and Green functions in the conformal basis. Suppressing again the w-dependence we have 


T(l'|l",2") 
r(2', 3"|4") 


M 3 ,M3.M3 >0 

M2,M 2 'M 2 >0 


r(5 // ,6"|l) = ]T r M' 1 M' 1 'lMl^M 1 ( 1 ) jB M' 1 ( 5// )^M' 1 '( 6,/ ), 

Mi,m'i-m" >° 

g(i"|6 ,/ ) = X] 

M4>0 

G'(2"|3") = Y Gm 5 ^ 5 (2 / ')E; 5 (3' / ), 

M5>0 

G(4"|5 // ) = Y G M 6 ^ 6 (4")^ 6 (5"). 

M6>0 

Integrations over the double primed coordinates will give 6 ^-functions in /j. So finally we 
are left with 6 summations over m. i=l,... 6 . The result can be presented in the form (with 
w-dependence suppressed) 


T(T,2'|1) = Y r M2,M3|M 1 ^M 1 (l)^(l , )^M3(2'), 

Ml,M2,M3>0 


(120) 
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where 


p _ p(0) 

M2?M3|Ml M2 5 M 3 |Ml 


8 gfjyg 

7T 2 


, \ ' 

2^? w ^ 


^6 |/M,M5 ^M2 ,M5 I, 


5 r M6,M4|Ati ^>4^>5 ^>6 


( 121 ) 

is the desired vertex part in the conformal basis. Its dependence on the c.m. coordinates Rj, 
i = 1,2,3, is determined according to Eq. (104). The bare vertex T^^^ is here given by 
the formula analogous to (107). 

In (121) summations over m, i = 4,5,6, include integrations over c.m. coordinates Rj, 
i = 4, 5,6. In the integrand the R dependence comes only from the vertex parts and in its 
turn is defined by Eq. (104). Thus the expression 


L hih2hz\h±hsh§ 


= R^Rrs 013 ^ 23 1 dR A dR b dR & R^R^R^R^R^R2 6 R^R^R^ • ( a.f .), ( 122 ) 

where all a’s are determined by formulas similar to (106) (see Appendix 2.), is independent 
of c.m. coordinates Ri, i = 1,2,3 and depends only on conformal weights. So the part of the 
vertex depending on conformal weights will satisfy an equation (with ui dependence restored 
in arguments) 


- - 8a 4 N 2 

r h2Mhl (u, u , ) = n(h 1 ,h2,h 3 )- s c 


7T 


7 \ ' 7 T 

‘2 7r j lfl 3 2_^ *'hG 1 hih 2 hzi\ / 1.4 h§ Hq 

h 4 ,h 5 ,h 6 >0 


^h 3 \h4,,h 6 ( UJ ^ u i)^h2te\h 6 { u, —u\,w—u')T h6 j li \hi( U} , ui)G hA {ui)Gh 5 (Lo'-ui)G h6 {u-ui). (123) 


7 Conclusions 

We have presented a formalism which allows to study interaction of pomerons in the QCD 
with N c —> 00 using the standard methods of quantum held theory. In particular we con¬ 
structed the Schwinger-Dyson equations which sum the so-called enhanced graphs and carry 
information of the 'physical’ pomeron as compared to the ’bare’ one. 

Conformal symmetry of the theory leads to certain simplifications. As a result we obtain 
a picture very similar to old Gribov local supercritical pomeron. The difference is reduced to 
an (infinite) number of pomerons with varying n = 0, ±2, ±4,... and more complicated form 
of the ’energy’ u: nv as a function of u, which plays the role of the pomeron momentum in the 
old theory, and of the bare triple pomeron vertex, which now depends both on n and u. If 
however one selects the supercritical pomeron with n = 0 and small values of u the formal 
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similarity becomes almost complete, since the energy becomes a quadratic function of v and 
the bare vertex then reduces to a well-know constant [17]. Unfortunately with this similarity 
also the problems of the old theory, mentioned in the Introduction, return together with the 
question of the internal consistency of the model. At present we do not know the answer to 
this question and leave it for future studies. 
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9 Appendix 1. Color factors 


Let the color wave function of a pair of gluons be | ab), where a,b = 1, ...Af? — 1. Then the 
vacuum color state is obviously 


where 


| 0 > = 




\aa) = P\ab) 


P = 




(124) 


(125) 


is the projector onto the vacuum color state. 

The color structure of the vertex for the transition 2—>4 reggeized gluons, with initial and 
final color variables a\, b\ —*► a-z, b- 2 , 03 , 63 is given by the expression 


-y _ ja\a2C jd>2d jda^ejeb^bi 


(126) 


We want the projection of this color vertex onto the 3 vacuum color states formed by the 
gluons with colors Gq&i, 0262 and 0363 . Applying three corresponding projectors Pi, P 2 and 
P 3 we obtain 


P2P3CP1 = 


(IV 2 - I ) 3 / 2 


j?aia,2Cj?ca2dfda^ejteasai _ 


JV 2 


(iv 2 - 1) 3 / 2 


^aid^aid 


N 2 




~ Nr 


(127) 

Note however that the quarks quark loop which represents the external source has its color 
factor 5 ab = s/N** - IP ~ N C P , so that each external source contributes a factor N c . 
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10 Appendix 2. Powers in Eq. (122) 

a 12 = ~ + ^(ni + n 3 - n 2 ) + i{v 2 -v$- vi) 

ai 3 = ~ + ^(ni + n 2 - n 3 ) + i (^ 3 - z / 2 - v\) 

«23 = - ^(ni + n 2 + n 3 ) + + z / 2 + u 3 ) 

a 34 = + ^(n 4 - n 5 - n 3 ) + i(v 5 - z / 4 + ^3) 

«35 = + 7^5 - n 4 - ra 3 ) + i(z/ 4 - z's + ^3) 

«45 = + ^(n 4 + n 5 + ra 3 ) - i(z / 5 + z / 4 + z/ 3 ) 

ai 6 = + ^(ni + n 4 - n 6 ) + i(z / 6 - z / 4 + z/ 4 ) 

ai 4 = + ^(ni + n 6 - n 4 ) + i(z / 4 - z / 6 + *h) 

a 46 = - ^(n 4 + n 4 + ra 6 ) + i(v 6 + ^4 + *h) 

a 2 6 = “ + ^(ra 2 + n 5 - ra 6 ) + i(^e - ^5 - ^2) 

a 25 = + ^(n 2 + n 6 - n 5 ) + i(z / 5 - i / 6 - z-'2) 

«56 = - ^(n 2 + n 5 + n 6 ) + i{v 6 + ^5 + ^2) 
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